Period of the 7-ray staggering in the 150 Gd superdeformed region 
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It has been previously proposed to explain 7-ray staggerings in the deexcitation of some superde- 
formed bands in the 150 Gd region in terms of a coupling between global rotation and intrinsic vortical 
modes. The observed 4h period for the phenomenon is suggested from our microscopic Routhian 
calculations using the Skyrme SkM* effective interaction. 
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This brief paper completes the theoretical investiga- 
tion of the coupling between global rotation and intrinsic 
vortical modes proposed in Ref. B as a tentative expla- 
nation for the rather rare band staggering observed in the 
decay of some superdeformed bands. Indeed while such a 
phenomenon was first claimed in 149 Gd [pj, then in 194 Hg 
H| , and possibly in the A ~ 130 superdeformed region j| , 
its existence has been confirmed and extended to a cou- 
ple of neighboring nuclei in the first case [0 and ruled 
out in the second case ||. Various theoretical explana- 
tions have been proposed |7Hllf| besides the one which 
we have discussed in Ref. [pT. When making an attempt 
to describe such data, one should address the three fol- 
lowing questions: (i) What is the mechanism at work? 
(ii) Why is this phenomenon so scarce and what makes 
it appear where it is observed (A ~ 150)? (iii) What is 
tuning the period of the staggering? 

While some answers have been provided in our previ- 
ous papers |I||l2|-|l4| to the two first questions, we aim 
here at addressing the third one. In [[l|, a staggering in 
transition energies within the yrast band was shown to 
appear in cases where the relevant collective energy is 
quadratic in two quantized quantities. A particular re- 
alization of the latter, well suited to the description of 
fastly rotating superdeformed states, corresponds to the 
parallel coupling of global rotation and intrinsic vorti- 
cal modes in ellipsoidally deformed bodies, known after 
Chandrasekhar jl5| as S ellipsoids. In this case, the two 
commuting operators are the projections on the quantifi- 
cation axis of the angular momentum operator and of the 
so-called Kelvin circulation operator (see, e.g., Ref. [fl6||), 
hereafter called / and J, respectively. Whereas it is triv- 
ial to show that the Kelvin circulation operator satisfies 
the usual commutation relations of an angular momen- 
tum, its consideration as a quantity which is approxi- 
mately a constant of the motion is a basic assumption of 
our collective model. Its exact amount of violation would 
deserve a specific microscopic study. A self-consistent de- 
scription of such a coupling can be made upon generaliz- 
ing the Routhian approach [T^ ] , amounting thus to solve 
the following variational problem: 

8{H — ClI — u)J) = 0, (1) 

where H is the microscopic Hamiltonian, Q (u>) an angu- 



lar velocity associated with the global rotation (intrinsic 
vortical) mode. This approach was first investigated in 
Ref. p2) within a simple oscillator mean field approxi- 
mation. There, the assumption of an energy which is 
quadratic in (f2,a;) or equivalently in (I, J) has been 
shown to be rather well satisfied. Recently, fully self- 
consistent solutions of the above variational problem have 
been made possible Jl8| upon using the standard SkM* 
force [^9). Details about such calculations and some of 
their results will be discussed below. 

In another paper |l3|| , a physical analogy stemming 
from the well-known similarity between the motion of 
a charge in a magnetic field and of a mass in a rotating 
frame has been established. It relates this staggering phe- 
nomenon with the observation of persistent currents in 
mesoscopic conductor or semiconductor rings as a man- 
ifestation of an inherent Aharonov-Bohm phase. Apart 
from its obvious physical appeal this consideration has 
set a framework in which one is able to understand the 
scarcity of both phenomena. It results from the necessity 
of securing a sufficiently low level of a specific damping 
which happens to yield in the staggering case a condition 
on the width of the superdeformed state in the relevant 
collective variable (e.g., the usual axial quadrupole defor- 
mation (3) . It was deduced from microscopic calculations 
of the associated mass parameters [B, using the usual 
D1S Gogny effective force pC] , that such a condition of 
existence was generally not met for Ce and Hg superde- 
formed states as well as for such states in Gd isotopes 
but for 150 Gd and possibly 148 Gd. 

As a consequence of all these studies the explanation 
for the staggering phenomenon suggested in Ref. is 
of course neither exempted from a priori questions nor 
deemed as being the only possible one, yet it is conforted 
as a rather likely candidate. However, one point remains 
to be clarified concerning the Ah period of the staggering. 
In Refs. Jl];[l2| this particular period was associated with a 
ratio of / and J values close to 2 for the considered states. 
From either semiclassical estimates of the relevant iner- 
tia parameters jlj or from actual microscopic calculations 
in the harmonic oscillator mean field approximation |l2j ] 
such a ratio is much more consistent with hyperdeforma- 
tion than with the actual superdeformation of the nuclear 
states. It is on this point that the new self-consistent ap- 
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proach of Ref. |18| brings some interesting insight. 

Here the generalized Routhian variational problem is 
solved within the Hartree-Fock approximation. Numer- 
ical codes breaking the time-reversal and axial symme- 
tries as requested by the considered physical problem, 
determine the single-particle wave functions cither at the 
nodes of a spatial mesh [2^] or as resulting from an ex- 
pansion on a suitably chosen truncated basis. In the lat- 
ter case, all approaches so far p2|-p4| , to the best of our 
knowledge, have used a triaxial basis. In our calculations 
an alternative method, shown to be less time consuming 
in most cases, has been developed where the expansion 
is made on an axial basis. Here the dependence of vari- 
ous fields and densities in terms of the angular variable 
of the cylindrical coordinate system is handled by conve- 
nient Fourier expansions. 



SkM* 
150 Gd 



80 - 



40 



50 



60 



/ (h) 



FIG. 1. Dynamical moment of inertia J^ 1 as a function of 
the angular momentum I. 

Self-consistent solutions of the usual Routhian prob- 
lem (i.e. without constraint on the Kelvin circulation, 
namely, for lu — 0) have been first performed for the 
150 Gd nucleus. As usual to get a solution corresponding 
to a given (quantized) value of /, one solves the problem 
iteratively so as to determine the angular velocity con- 
straint Q able to reach the requested angular momentum. 
For such solutions, one can estimate the dynamical mo- 
ment of inertia by numerical differentiation. As suggested 
in Ref. |2J|, among the possible expressions for j( 2 \ we 
choose for numerical stability reasons to evaluate 
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Such derivatives have been determined from a two-point 
formula involving thus for each studied value of the angu- 
lar momentum three different Hartree-Fock calculations 
(differing typically around I = 50H by A?if2 ~ 0.07 keV). 
As discussed, e.g., in Ref. |jl2), even in the absence of any 
constraint on the intrinsic vortical currents, the Kelvin 
circulation operator takes a finite value which is easily 
estimated from our solutions. The vector components of 



this operator are defined (see, e.g., Ref. |ll|) by its action 
on single-particle wave functions as 



h ^ c 7 d 
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where e a /3 7 is the completely antisymmetrical third-rank 
tensor, x a the a component of the particle position vec- 
tor, and c a the corresponding length scale factor. These 
factors are deduced from the values of the quadrupole 
tensor calculated from our variational solutions upon 
making an ellipsoidal shape approximation. 

Some results of our Routhian calculations (uj — 0) for 
a range of / values from / = 407i to / = 66h, including 
states of relevance for the observed superdeformation, are 
displayed in Fig. |l|. It shows the variation of the moment 
of inertia as a function of /. As a matter of fact, it 
fits rather well, as it should, with the pure Hartree-Fock 
part of the results obtained in Ref. [|26f with the same 
interaction. The calculated yrast value J yra st(-0 of the 
Kelvin circulation as a function of / falls very nicely on 
the straight line: 



•WCO ^0.8J+1.0?L 
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This is not at all surprising insofar as the quadratic ap- 
proximation for the collective energy discussed in Ref. 
p"2| is valid over the whole considered range of values of 
/. However, this relation makes it very clear that around 
/ = 50h, e.g., the ratio I/J is indeed very far from the 
value of 2 which would lead to the observed staggering 
period. 
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FIG. 2. Total energy £i{J) as a function of the Kelvin cir- 
culation J for angular momentum values around I = 50h. 

As a result of the continuity equation, tangential in- 
trinsic vortical excitations yield phase-space modifica- 
tions amounting only to a momentum redistribution [fl2[ . 
In that respect they are indeed quite similar to pairing 
correlations. This has been, for instance, illustrated from 
another point of view in Ref. p7| . There, current pat- 
terns as functions of a fixed pairing gap display indeed 
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the same type of variation as classical S'-ellipsoid velocity 
fields with respect to the angular velocity of the intrinsic 
vortical modes. Therefore one may consider the latter 
modes as a collective model translation of pairing corre- 
lations in a somewhat similar fashion as small amplitude 
vibrational collective modes can model random phase ap- 
proximations (RPA) correlations. 

In order to implement this type of excitation on top of 
Hartree-Fock rotational solutions we have solved the gen- 
eralized Routhian problem with two constraints (ft and 
uj, both ^ 0). Indeed we have made calculations for fixed 
values of / upon varying J. Clearly the single constraint 
corresponds to the yrast state (up to the quantization of 
J of course). It is therefore no surprise to find it at the 
minimum of a somewhat parabolic energy curve £i{J) as 
demonstrated in Fig. |[ Note in passing that the general 
pattern of such energy curves £i(J) for various values of 
I is consistent with a quadratic dependence of the total 
energy in / and J as assumed in Rcf. [0 and calculated 
in the simple model case of Ref . [n~2| . 
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FIG. 3. Variations of the two relevant angular velocities Q. 
and ui as functions of the Kelvin circulation J for I = 50ft 
solutions. 

Now to make excursions, for a given value of /, out of 
the yrast solution, one has to perform a two Lagrange 
multipliers search, whose result is examplified in Fig. y 
for the I = 50h case. It is rather significant that to get 
J values smaller than the yrast value, one should add 
a counterotating intrinsic vortical mode. This can be 
explained by using the quadratic approximation for the 
total energy. One finds 



I = CVL + Buo 



J = Btt + Aw 



(5a) 
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where A, B, and C are (positive) inertia parameters de- 
fined in Ref. (L2|. From the above, one finds trivially 



where the coefficient of u> is found to be positive for well- 
deformed nuclei as easily seen from the semiclassical es- 
timates of Refs. fUJl^. Therefore one can check that 
starting from fl > for the yrast state, in order to de- 
crease J, one diminishes (from its vanishing value) the u> 
velocity while increasing £1 (keeping / constant) so that 
one gets wft < 0. Conversely the same reasoning yields 
uiQ > when increasing J away from its yrast value. 

Pairing correlations act against the global rotations 
(see, e.g., Ref. |27|]). In our collective model account of 
such correlations, this corresponds to a product of an- 
gular velocities such that wO < 0. Consequently start- 
ing from a noncorrelated (Hartree-Fock) solution, the in- 
clusion of pairing correlations will tend to decrease the 
Kelvin circulation value (see Fig. ||). In Fig. |J, the 
dynamical moments of inertia are plotted for three 
values of / around 50?i as functions of J. One sees that 
pairing correlations will indeed increase from its 

Hartree-Fock value as actually found in the Hartree-Fock- 
Bogoliubov calculations of Ref. p^j. In this paper, the 
authors have shown upon using simple yet realistic pair- 
ing matrix elements that the correlations raise the mo- 
ment of inertia to the vicinity of the experimental value 
(typically ~ 90?i 2 MeV -1 ). It is striking that when 
constraining the intrinsic vortical mode to obtain this 
value for J^ 2 \ one gets a Kelvin circulation of ~ 25?i 
which fulfills precisely the I / J ratio condition of being 
close to 2. The latter provides a Ah period for the oscil- 
lating behavior of the 7 transition energies in the 150 Gd 
region. 



5» 









SkM* 


- \ 






150 Gd 








/ = 48fi 








/ / = 50fi 








/ / I = 52fi. 



















J c 



AC~B 2 

i 

C 



(6) 



J (Ti) 

FIG. 4. Dynamical moment of inertia as a function of 
the Kelvin circulation J for angular momentum values around 
/ = 50/2. Dots correspond to yrast solutions. 

To confirm the above conclusion, it would be, of 
course, very interesting to perform variational gener- 
alized Routhian calculations within the Hartree-Fock- 
Bogoliubov approximation. We are currently working on 
it. Nevertheless, it seems to us very likely that our indi- 
rect estimate already leads to the conclusion that pairing 
correlations play the major role in fine-tuning the stag- 
gering period to what is experimentally observed. 
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